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ABSTRACT 

We use Type la supernovae studied by the High-Z Supernova Search Team to 
constrain the properties of an energy component which may have contributed to 
accelerating the cosmic expansion. We find that for a fiat geometry the equation 
of state parameter for the unknown component, ttx = -Px/Px? must be less than 
—0.55 (95% confidence) for any value of fl^a and is further limited to ttx < —0.60 
(95%) if Qui is assumed to be greater than 0.1 . These values are inconsistent 
with the unknown component being topological defects such as domain walls, 
strings, or textures. The supernova data are consistent with a cosmological 
constant {a^ = — 1) or a scalar field which has had, on average, an equation of 
state parameter similar to the cosmological constant value of —1 over the redshift 
range of 2; ~ 1 to the present. Supernova and cosmic microwave background 
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observations give complementary constraints on the densities of matter and the 
unknown component. If only matter and vacTium energy are considered, then the 
current combined data sets provide direct evidence for a spatially flat Universe 
with Qtot = flra + ^A = 0.94 ± 0.26 (la). 

Subject headings: supernovae — cosmology: observations and cosmic microwave 
background 



1. Introduction 

Matter that clusters on the scale of galaxies or galaxy clusters is insufficient to close the 
Universe, with conventional values near Qm = 0.2 ±0.1 (Gott et al. 1974; Carlberg et al. 
1996; Lin et al. 1996; Bahcall, Fan, & Cen 1997). Observations of distant supernovae provide 
credible evidence that the deceleration rate of the Universal expansion is small, implying 
that the total matter density, clustered or smooth, is insufficient to create a flat geometry 
(Garnavich et al. 1998; Perlmutter et al. 1998). Either the Universe has an open geometry or, 
if flat, other forms of energy are more important than matter. 

Large samples of supernovae analyzed by the High-Z Supernova Search collaboration 
(Riess et al. 1998a, hereinafter [Riess98]) and the Supernova Cosmology Project (Kim 1998) 
now suggest that the Universe may well be accelerating. Matter alone cannot accelerate the 
expansion; so if taken at face value, the observations demand an additional energy component 
for the Universe. While the vigorous pursuit of possible systematic effects (e.g. Hoflich, 
Wheeler, & Thielemann 1998) will be important in understanding these observations, it is 
instructive to see what they imply about the energy content of the Universe. 

The cosmological constant was revived to fill the gap between the observed mass density 
and the theoretical preference for a fiat Universe (Turner, Steigman, & Krauss 1984; Peebles 
1984) and also to alleviate the embarrassment of a young Universe with older stars (Carroll, 
Press, & Turner 1992). The cosmological constant is a negative pressure component arising 
from non-zero vacuum energy (Weinberg 1989). It would be extraordinarily difficult to detect 
on a small scale, but Qa = 1 — f2m could make up the difference between the matter density 
Jim and a flat geometry and might be detected by measurements on a cosmological scale. 
There are few independent observational constraints on the cosmological constant, but Falco, 
Kochanek, &: Muiioz (1998) estimated that Q,\ < 0.7 (95% confidence) from the current 
statistics of strong gravitational lenses. If the matter density is less than fim ~ 0.3, this 
limit is close to preventing the cosmological constant from making a fiat geometry. Further, 
a cosmological constant which just happens to be of the same order as the matter content at 
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the present epoch raises the issue of "fine tuning" (Coles & EUis 1997). A number of exotic 
forms of matter which might contribute to cosmic acceleration are physically possible and 
viable alternatives to the cosmological constant (Prieman & Waga 1998; Caldwell, Dave, & 
Steinhardt 1998). The range of possibilities can be narrowed using supernovae because the 
luminosity distance not only depends on the present densities of the various energy components 
but also on their equations of state while the photons we see were in flight. Here, with some 
simplifying assumptions, we consider the constraints that recent supernova observations place 
on the properties of an energy component accelerating the cosmic expansion. 



2. Observations 

The type la supernovae (SNIa) have been analyzed by the High-Z Supernova Search Team 
and described by Riess98, Garnavich et al. (1998), Schmidt et al. (1998), and Riess et al. 
(1998b). The full sample from Riess98 consists of 50 SNIa. Of these, 34 are dX z < 0.2 while 
the remaining 16 cover a range in redshift of 0.3 < z < 1.0. Six of the high-redshift events 
were analyzed using the "snapshot" method developed by Riess et al. (1998b). This innovative 
technique uses high-quality spectra to deduce information unavailable due to a poorly sampled 
light curve. While the errors estimated from the snapshot method are larger than those from 
direct light curve fitting, the snapshot sample provides a significant, independent set of SNIa 
distances. 

As shown by Philhps (1993), the fight curve decline rate of SNIa is correlated with the 
luminosity at maximum brightness of these exploding white dwarfs. This correlation has been 
cafibrated by Hamuy et al. (1996, the Ami5(S) method) and by Riess, Press, & Kirshner 
(1995, 1996, the Multi-Color Light Curve Shape or MLCS method which includes a correction 
for extinction), and both show that applying this correction to the SNIa Hubble diagram 
significantly reduces the scatter. Phillips et al. (1998) extended the ^mi^{B) approach to 
include an estimate of the extinction. In Ricss98, an improved version of the MLCS method 
is presented. Here, as in Riess98, we apply both MLCS and Am 15(5) (with extinction 
correction) techniques to the analysis to gauge the systematic errors introduced by different 
light curve fitting methods. 



3. Analysis 

The apparent brightness of a SNIa corrected for light curve decline rate and extinction 
provides an estimate of the luminosity distance, D^, from the K-corrcctcd observed magnitude, 
m = M + blogDi + 25, and the absolute magnitude, M, of SNIa. As described by Schmidt et 
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al. (1998) and Carroll, Press, & Turner (1992), the luminosity distance depends on the content 
and geometry of the Universe in a Priedmann-Robertson- Walker cosmology 

= "^^tA strmlM^l /"[E^^(l + ^0'^'+"^^ + ^k{l + z'fr'^'dA (1) 

Ho^J\ni:\ r -^0 Y J 

' sinh{x), if flk > 0; 
sinn{x) — I if = 0; (2) 

, sin{x), if Qjfc < 0, 

where fij are the normalized densities of the various energy components of the Universe 
and = ^ — Y^i^i describes the effects of curvature. The exponent n = 3(1 + a) defines 
the way each component density varies as the Universe expands, p oc a~", where a is the 
cosmic scale factor. For example, n has the value 3 for normal matter since the mass density 
declines proportionally to the volume. Alternatively, ccj is the equation of state parameter 
for component i defined as the ratio of the pressure to the energy density, ai — Pi/ Pi 
(sometimes denoted in the literature as w). The relation n = 3(1 + a) is easily derived from 
the conservation of energy equation in comoving coordinates (e.g. Weinberg (1972) equation 
15.1.21). In the most general case, the equation of state can vary with time in ways other 
than assumed here (as the sum of power laws in 1 + 2;), but we are limited by the quality and 
range of the supernova observations to consider only its average effect between the present 
and z < 1. The present-day value of the Hubble constant (Hq) and the absolute magnitude of 
SNIa (M) are primarily set by the low-redshift sample, which allows the high-redshift events 
to constrain the cosmological effects. This means that conclusions derived from SNIa are 
independent of the absolute distance scale. 

Gravitational lensing by matter distributed between the observer and supernovae at high 
redshift can affect the observed brightness of SNIa and induce errors in the estimate of their 
luminosity distances (Kantowski, Vaughan, & Branch 1995). For reahstic models of the matter 
distribution and Jim < 0.5, the most likely effect of the lensing is to make the supernovae at 
z — 0.5 about 2% fainter than they would appear if the matter were distributed uniformly 
(filled-beam) as shown by Wambsganss et al. (1997). Holz & Wald (1997) have shown that the 
magnitude of the effect also depends on whether the matter is distributed smoothly on galaxy 
scales or is clumped in MACHOS, but the error induced remains small when Q^i < 0.5. For 
simplicity, our calculations consider only the filled-beam case, however, the effect of assuming 
the extreme case of an empty-beam is shown by Holz (1998). 

There are a few known, and possibly some unknown, energy components that affect Dl. 
Ordinary gravitating matter, certainly has had some effect on the Universal expansion 
between z ^ 1 and now. Since the matter density scales inversely with the volume, am = 0, 
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and matter (baryons, neutrinos, and dark matter; formerly Earth, Air, and Water) contributes 
no pressure. Radiation (Fire in an earlier lexicon) (aj. = +1/3) dominated during a period in 
the early Universe but is negligible for z < 1. Equation 1 shows that for non-flat models the 
curvature term, Qk, contributes to the luminosity distance like a component with ak — —1/3, 
but additional geometrical effects as prescribed by equation 2 are also important. 

Other more speculative components have been proposed. A non-zero vacuum energy, 
is a popular possibility explored by Riess98 for this data set. Because the vacuum energy 
density remains constant as the Universe expands (that is, Pa oc a°), we have a a = —1. 
Topological defects created in the early Universe could also leave remnants that might 
contribute to the energy now. Networks of cosmic strings may be a natural consequence of 
phase transitions in the young Universe and if they did not intercommute would have an 
average effective ctg = —1/3 (Vilenkin 1984; Spergel & Pen 1997). A network of comoving 
domain walls would have an average equation of state parameter of —2/3 (Vilenkin 1985) 
while a globally wound texture would produce an at — —1/3 (Davis 1987; Kamionkowski & 
Toumbas 1996). 

Evolving cosmic scalar fields with suitable potentials could produce a variety of exotic 
equations of state with significant densities at the present epoch (Peebles & Ratra 1988; 
Prieman et al. 1995; Prieman & Waga 1998). Scalar fields could also produce variable mass 
particles (VAMPS) which would redshift more slowly than ordinary matter creating an 
effective ckvamp < (Anderson & Carroll 1998). These fields may evolve over time and would 
produce an interesting variety of cosmic histories. Our goal is modest: we only hope to 
constrain the average a over the range where SNIa are presently observed. 

To simplify the analysis, we assume that only one component affects the cosmic expansion 
in addition to gravitating matter. Because the origin of the acceleration is unknown, we will 
refer to this as the "X" component with a density of fl^ and equation of state of — Oi^Px- 
Caldwell, Dave, & Steinhardt (1998) have dubbed the unknown component "quintessence" as 
the other four essences have ahcady been employed above. We assume that the Universe on 
very large scales is accurately described by general relativity and that the "X" component 
obeys the null energy condition (NEC). The NEC states that, for any null vector v'^, the 
energy-momentum tensor satisfies T^^v^^v" > (see, e.g., Wald 1984). This is the weakest of 
all conventional energy conditions, and should be satisfied by any classical source of energy 
and momentum including those discussed above. In a Robertson- Walker metric, the NEC is 
equivalent to requiring Px + -Px > 0. The NEC therefore restricts the energy density of the 
unknown component to be positive for > —1 and negative when ccx < — 1 while the energy 
density of the cosmological constant (ccx = —1) is unconstrained. 
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Fig. 1. — The joint probability distributions for fim and the density of the unknown component, 
fl^, based on the SNIa magnitudes reduced with the MLCS method. Four representative values 
of the equation of state parameter, ctx, are shown. See Riess98 for the distribution when 
ax = -1- 



4. Results 

First, we fix the equation of state of the unknown component and estimate the probability 
density function for the parameters Q^, ^m, and Hq given the observed SNla distance moduli. 
The joint likehhood distributions are then calculated in the same way as by Riess98 and 
shown for representative values of ccx in Figure 1. Here we integrate over all possible Hq with 
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the prior assumption that all values are equally likely. For ccx < —0.7, the derived constraints 
are similar to those found by Riess98 for a cosmological constant (ctx = — 1). However, as the 
equation of state parameter increases, the major axis of the uncertainty ellipses rotates about 
a point on the f2x = line. For an accelerating Universe, the pivot point is on the negative 
f2m side. When cix > —0.4, the "X" component could not reproduce the observed acceleration 
and all of these models give a poor fit to the observed SNIa data: the best fit occurs for a 
completely empty Universe. 

Next, we allow the equation of state to vary freely but restrict the densities to r^m + rix < 1, 
or open models. We then integrate over all possible values of assuming a uniform prior 
distribution to provide the joint probability for and the matter density. From the NEC 
we must include regions where ctx < —1 and is negative, but these are unable to produce 
an accelerating Universe so they have a very low probability. For open models, highest 
joint probabilities are confined to a region bounded by —1.0 < ccx < —0.4 and < 0.2. 
If we consider any value of 0,^ equally hkely, then ctx < —0.47 for the MLCS method and 
CKx < —0.64 for the Ami5{B) results with 95% confidence. 

Finally, we consider fiat models for the Universe. The joint probability between the 
equation of state parameter and the matter density for fl^ + Qx = 1 is shown in Figure 2. The 
two cases are for the MLCS and the Ami5{B) hght curve fits and demonstrate that the two 
methods for deriving luminosity from light curves provide consistent constraints. Note that in 
the flat case, the NEC allows cKx < —1 only when fl^ > 1 which has an insignificant probability 
and is not plotted. These plots can be compared to pioneering calculations by Turner & 
White (1997) and White (1998) which used smaller supernova samples. The improved SNIa 
data favor acceleration and support both a low Q^a and a small value of a^- Integrating the 
probability over all values of flra assuming a uniform prior shows that < —0.55 for MLCS 
and ccx < —0.63 for Ami5{B) (95% confidence). If we assume fl^ > 0.1 then the fimits tighten 
to CKx < —0.60 (MLCS) and ccx < —0.69 (Ami5(i?)). For matter densities near ~ 0.2 favored 
by galaxy cluster dispersions, the most probable equation of state parameters are between 
—0.7 and —1.0. These results disfavor topological defect models such as domain walls (90% 
confidence) and eliminate strings and textures (99% level) as the principal component of the 
unknown energy. The cosmological constant, or a form of quintessence that resembles it for 
z < 1, is supported by the data. Constraints that refer to higher redshift are needed to narrow 
the range of possible models. 
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5. Other Constraints 

High-Z SNIa observations combined with the cosmic microwave background (CMB) 
anisotropy angular power spectrum provide complementary constraints on the densities of 
matter and the "X" component (White 1998; Tegmark et al. 1998). Details of the CMB power 
spectrum depend on a large number of variables, but the angular scale of the first acoustic 
peak depends primarily on the physics of recombination and the angular diameter distance to 
the surface of last scattering (White & Scott 1996; White 1998). Rather than fit the power 
spectrum in detail, we have restricted our attention to the location of the first acoustic peak 
as estimated from current CMB experiments (Hancock et al. 1998). This is a rapidly moving 
experimental field, and new results will surely supersede these, but they illustrate the power 
of combining the supernova data with the CMB. We employ the analytic approximations of 
White (1998) to determine the wavenumber of the acoustic peak at recombination, and those 
of Hu & Sugiyama (1996) to determine the recombination redshift; thus we assume adiabatic 
fluctuations generate the anisotropy. In addition, we have ignored reionization and flxed the 
number of neutrino species at three, as well as assuming only scalar modes, with a spectral 
index n — 1. K thorough treatment of this problem would allow all of these parameters to vary 
and integrate the probability over all possible values (that is, marginalize over them); however, 
this would be very time-consuming, even with the fast CMB code of Seljak & Zaldarriaga 
(1996), and disproportionate to the precision of the current data. A large exploration of the 
parameter space involved (though lacking a full variation of Vt^ can be found in Bartlett et 
al. (1998) and Lineweaver (1998). 

Our calculation determines the angular scale multipole of the flrst acoustic peak for a 
grid in a three-dimensional parameter space of [VLmi^KiHo), where we explicitly allow for 
open, flat, and closed universes with and without a cosmological constant. We also employ 
the additional constraint on the baryon density Vtfjh^ = 0.024 {h = i/o/100 kms^^Mpc"^) 
derived from the primordial deuterium abundance and nucleosynthesis (Tytler, Fan, & Buries 
1996). Other estimates of the baryon fraction (see Fugikita, Hogan, & Peebles 1998) could be 
used, but the location of the peak depends only weakly on this parameter. Where possible, 
we checked these calculations with numerical integrations (Seljak & Zaldarriaga 1996) and 
conflrmed that the peak locations agree to ^ 10%, which is adequate for this exploration. 

Following White (1998), we combine the predicted peak location with the observations 
using a phenomenological model for the peak (Scott, Silk, & White 1995). Recent CMB 
measurements analyzed by Hancock et al. (1998) give the conditional hkelihood of the flrst 
acoustic peak position as Zpeak = 263l!l94^, based on best-flt values of the peak amplitude 
and low multipole normalization. Rocha et al. (1998) have provided us with a probability 
distribution function for the flrst peak position based on marginalizations over the amplitude 
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and normalization which is a more general approach than by Hancock et al. The Rocha et 
al. function gives /peak = 284184^ which is only a small shift from the value derived using the 
conditional likelihood method. Wc then marginalize the likelihood in our three-dimensional 
parameter space over Hq with a Gaussian prior based on our own SNIa result including our 
estimate of the systematic error from the Cepheid distance scale, i^o = 65 ± 7 kms^^ Mpc~^ 
(Riess98). It is important to note that the SNIa constraints on (f2in,^A) are independent of 
the distance scale, but the CMB constraints are not. We then combine marginalized hkelihood 
functions of the CMB and SNe la data. The result is shown in Figure 3. Again, we must 
caution that systematic errors in either the SNIa data (Riess98) or the CMB could affect this 
result. 

Nevertheless, it is heartening to see that the combined constraint favors a location in this 
parameter space which has not been ruled out by other observations, though there may be 
mild conflict with constraints on from gravitational lensing (Falco, Kochanek, & Muiioz 
1998). In fact, the region selected by the SNIa and CMB observations is in concordance with 
inflation, large-scale structure measurements, and the ages of stars (Ostriker & Steinhardt 
1995; Krauss & Turner 1995). The combined constraint removes much of the high Qm, high 
f^A region which was not ruled out by the SNIa data alone, as well as much of the high Q^, 
low f^A region allowed by the CMB data alone. The combined constraint is consistent with 
a flat universe, as Qtot = ^Im + = 0.94 ± 0.26 for MLCS and 1.00 ± 0.22 for Amis (5) 
(la errors). The enormous redshift difference between the CMB and the SNIa makes it 
dangerous to generahze this result beyond a cosmological constant model because of the 
possible time-dependence on ccx- But for an equation of state flxed after recombination, the 
combined constraints continue to be consistent with a flat geometry as long as ccx ^ —0.6. 
With better estimates of the systematic errors in the SNIa data and new measurements of the 
CMB anisotropy, these preliminary indications should quickly turn into very strong constraints 
(Tegmark et al. 1998). 
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Fig. 2. — The joint probability distributions from SNIa for and the equation of state 
parameter, a^, assuming a flat spatial geometry (fim + Qx = !)• The top panel uses supernova 
distances from the MLCS method combined with supernovae reduced using the snapshot 
method, while the bottom panel is from the Amis (5) technique plus snapshot results. The 
vertical broken line marks the matter density estimated from galaxy cluster dynamics. 
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6. Conclusions 

The current results from the High-Z Supernova Search Team suggest that there is 
an additional energy component sharing the Universe with gravitating matter. For a flat 
geometry, the ratio of the pressure of the unknown energy to its density is probably more 
negative than —0.6. This effectively rules out topological defects such as strings and textures 
as the additional component and disfavors domain walls as that component. Open models are 
less constrained, but favor ccx < —0.5. Although there are many intriguing candidates for the 
"X" component, the current SNIa observations imply that a vacuum energy or a scalar field 
that resembles the cosmological constant is the most likely culprit. 

Combining the SNIa probability distribution with today's constraints from the position 
of the first acoustic peak in the CMB power spectrum provides a simultaneous observational 
measurement of the densities of matter and of the unknown component. Using CMB data 
from Hancock et al. (1998) and following the analysis by White (1998) the result favors a flat 
Universe, with fltot = 0.94 ± 0.26, dominated by the "X" component for fa —1. Given the 
rapid improvement in both the study of SN la and the CMB, we can expect more powerful 
inferences about the contents of the Universe to follow. 
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Fig. 3. — The combined constraints from SNIa and the position of the first Doppler peak of the 
CMB angular power spectrum. The equation of state parameter for the unknown component is 
CKx = —1, hke that for a cosmological constant. The contours mark the 68%, 95.4%, and 99.7% 
enclosed probability regions. 
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